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ABSTRACT 

The transient compressible one-dimensional vapor 
flow dynamics in a heat pipe is modeled. The numerical 
results are obtained by using the implicit non-iterative 
Beam-Warming finite difference method. The model is 
tested for simulated heat pipe vapor flow and actual vapor 
flow in cylindrical heat pipes. A good comparison of the 
present transient results for the simulated heat pipe vapor 
flow with the previous results of a two-dimensional 
numerical model is achieved and the steady state results are 
in agreement with the existing experimental data. The 
transient behavior of the vapor flow under subsonic, sonic, 
and supersonic speeds and high mass flow rates are 
successfully predicted. The one-dimensional model also 
describes the vapor flow dynamics in cylindrical heat pipes 
at high temperatures. 

NOMENCLATURE 

0^ specific heat at constant pressure, J/kg-K 
C y specific heat at constant volume, J/kg-K 
D y diameter of the vapor space, m 
E t total energy of the vapor per unit volume, p (C y T + 

1/2 U 2 ) 

2 

f friction coeflicient at the wall, r Q /pU 

f g friction coefficient at the exit of the evaporator 

f^ friction coefficient for the fully developed turbulent 
flow 

h enthalpy, J/kg 

IIj£ latent heat of vaporization, J/kg 

K thermal conductivity of the vapor, W/m-K 

L length of the heat pipe, m 

L a length of the adiabatic section, m 

L length of the condenser, m 

L e length of the evaporator, m 

M Mach number, U//fRT 


n time step 

o 

P vapor pressure, N/ra 

P c reference pressure for the Clausius-Clapeyron 
2 

equation, N/m 

O 

P fi pressure at the evaporator end cap, N/m 

2 

P g k sink pressure, N/m 

2 

P g0 source pressure, N/m 

Q heat input, W 

o 

q heat flux, W/m 

R gas constant, J/kg-K 

Re axial Reynolds number, pUD y //x 

Re Q radial Reynolds number at the wall, pV Q D y //x 

t time, s 

T vapor temperature, K 

T c reference temperature for the Clausius-Clapeyron 
equation, K 

U axial 'velocity, m/s 

V Q radial velocity at the wall, m/s 

x coordinate in the axial direction 

x t . initial location of the transition region, m 

x dimensionless location of the transition region, (x - 

* tii )/n 

Greek Symbols 

7 ratio of specific heats, C p /C y 

£t time increment, s 

fa distance between nodes, m 

T fraction of transition to turbulent flow 

2 

/x dynamic viscosity of the vapor, N— s/m 

2 

p density of the vapor, kg/m 

t shear stress 

n x r = 3/4 ~ x r = 1/4 
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o properties of the injected or extracted fluid at the 

wall 

x derivative with respect to x 
INTRODUCTION 

In a heat pipe, the change of phase of the working 
fluid in the dosed system is used instead of a large 
temperature gradient to transport a large amount of energy. 
The attention of many sdentists has focused on this unique 
phenomenon since the concept was introduced. The vapor 
from the evaporator carries energy to the condenser, so the 
vapor flow in the core region of the heat pipe plays an 
important role in transferring energy from source to sink. 

Many researchers have studied the steady one- 
dimemional compressible (Levy, 1968; Brovalsky et aL, 
1976; Faghri, 1988; and Jang, 1988) and the steady 
two-dimensional vapor flow in heat pipes (Bankston and 
Smith, 1972; Tien and Rohani, 1974; Ooijen and 
Hoogendoorn, 1979; Faghri, 1986; and Faghri and Parvani, 
1988). The common cross sections of the vapor space are 
circular, rectangular (Jang, 1988; Ooijen and Hoogendoorn, 
1979), and annular (Faghri, 1986; Faghri, 1988; and Faghri 
and Parvani, 1988) and are chosen based on the particular 
application. The heat flux distributions on the surface of 
the evaporator and condenser are uniform except for those 
presented by Jang (1988). Recently, the transient 
two-dimensional compressible simulated vapor flow in heat 
pipes (Bowman, 1987; and Bowman and Hitchcock, 1988) 
was solved numerically, and the experimental data for the 
steady state simulated heat pipe vapor flow was obtained by 
Bowman (1987). However, the vapor flow was not actual 
vapor flow in a heat pipe and the numerical and 
experimental data were presented only in terms of the 
pressure. The transient two-dimensional compressible 
vapor flow in a heat pipe with a rectangular cross section 
also was studied numerically by Issacci et al. (1988), but 
only the axial and radial velocity profiles were presented. 
No comparison with the existing experimental data was 
made and the transient behavior of vapor flow in a heat pipe 
was not described. 

During the start-up of high temperature heat pipes, 
the extremely small density of the vapor causes the vapor 
flow to attain sonic and supersonic velocities for a relatively 
small heat input. Thus, the correct description of the 
transient vapor flow is essential to predict the successful 
start-up and to estimate the overall performance of the 
entire heat pipe. The governing equations for the vapor 
flow as well as those for the wall and wick regions should be 
solved simultaneously. Also, the development of the 
one-dimensional transient model for the vapor flow has 
been suggested due to the large amount of computer time 
required for the two— dimensional model (Bowman, 1987; 
and Bowman and Hitchcock, 1988). For this purpose, a 
simple and efficient mathematical model is desired for each 
region. All previous one— dimensional models for vapor flow 
are for the steady state condition and the viscous dissipation 
was neglected, which is important for high temperature 
applications. 

This paper describes the mathematical model and the 
numerical method of solution for the transient compressible 
one— dimensional vapor flow in the heat pipe. A comparison 
of the numerical results with the simulated transient 
two-dimensional numerical results and experimental data 
for the steady state given by Bowman (1987) is also 
presented. In addition, the numerical results from the 
present model for the actual vapor flow in the cylindrical 
high temperature heat pipe are compared with the 
experimental data obtained at the steady state by 
Ivanovskii et al. (1982). 


MATHEMATICAL MODELING 


The one-dimensional transient compressible vapor 
flow is considered to predict the vapor flow dynamics in the 
heat pipe. Even though a uniform velocity is used, the 
friction at the interface is incorporated by using the 
expressions for the friction coefficients which are found from 
the two-dimensional numerical results given by Bowman 
(1987). The viscous dissipation in the vapor region is 
included and the vapor is assumed to be a perfect gas. The 
governing equations for the vapor flow in the heat pipe with 
negligible body forces are formulated by using the principles 
of the conservation of mass, momentum, and energy in a 

2 

control volume of cross-sectional area, x D y / 4 , and width, 
dx. 

The governing equations are written in a compact 
vectorial form as follows: 


®+^-k f ‘ d ' d .) +g 
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(4) 




(5) 



where V Q (x) refers to the velocity at the wall with a 

positive value for injection and a negative value for suction. 

The equation of state is employed to relate the 
density, pressure, and temperature in the vapor space as 
follows: 

P = pRT (6) 


For the simulated heat pipe vapor flow, the temperature 
was evaluated by using the equation of state because a 
change of phase was not involved. For the actual vapor flow 
in the cylindrical heat pipe, the Clausius-Clapeyron 


2 



ORIGINAL PAGE !S 
OF POOR QUALITY 


relationship was used to predict the saturation temperature 
of the vapor from the pressure as given by 


T = 


1 K' ]n P 

T c - *7^ *1 


(7) 


The known boundary conditions at the ends of the 
heat pipe are as follows: 


and output is important. To eliminate this difficulty, terms 
in braces {1 in equation (5) are replaced by using the heat 
flux appliea on the surfaces of the evaporator and condenser 
as follows: 

r v 2 (x) 

<1 = W*) [ h 0 W + — 5 — ] (13) 

FRICTION COEFFICIENTS 


U = 0 at 

x = 0 and L 

(8) 

S = 0 at 

x = 0 and L 

(9) 


The conditions for the density and pressure at the ends of 
the heat pipe are unknown, so physically realistic boundary 
conditions should be derived. In general, near the ends of 
the heat pipe the mass flow rate is small, so the axial 
gradients of the pressure and density are small. In the region 
adjacent to the exit of the evaporator, the variations of the 
pressure and density are large. Thus, the boundary 
conditions for the pressure and aensity can be assumed as 
follows: 


d 

o 

II 

x — 0 and L 

(10) 

d 

O 

ii 

x = 0 and L 

(11) 


Bowman (1987) introduced a correlation between the 
mass flux at the wall and the pressure drop across a porous 
tube wall based on experimental measurements for the 
simulated heat pipe vapor flow. The mass flux (p Q V 0 ) at 

the wall for the simulated heat pipe vapor flow was 
evaluated by using this correlation: 

A(P 2 * * * * * ) = 3.639 x 10V o V o ) 2 + 1.7015 x 10 8 (/, o V o ) (12) 

2 

where A(P ) is the absolute value of the difference between 

the square of the uniform source pressure and the square of 

the vapor pressure in the blowing section. In the suction 

o 

section, A(P ) is the absolute value of the difference 
between the square of the vapor pressure and the square of 
the uniform sink pressure. A change of phase of the working 
substance was not involved. The uniform source 
temperature of 300 K was used for his experiment, but the 
sink temperature was not specified. To evaluate the terms 
in braces {} in equation (5), the source temperature is used 
at the blowing section and the vapor temperature is 
employed at the suction section. 

Since the working fluid changes phase at the 
vapor-liquid interface in actual heat pipes, the temperature 
at the interface is the saturation temperature, but the 
temperature in the vapor space may be quite different from 
the saturation temperature for high temperature heat pipes. 
For the one-dimensional model, the properties are the 
area-averaged properties so that the temperature in the 
vapor space is not the interface temperature but is also not 
the saturation temperature. The vapor temperature can be 
evaluated from the energy equation and the saturation 
temperature corresponding the vapor pressure can be 
obtained from the Clausius— Clapeyron relationship. 
However, this saturation temperature is not the actual 
interface temperature either. The correct estimation of the 
terms in braces {} in equation (5) is uncertain due to using 
the area- averaged properties. Since a heat pipe is a closed 
system, the application of the correct values of heat input 


Bowman (1987) measured the turbulence intensity in 
the simulated heat pipe to determine the characteristics of 
the vapor flow. For axial Reynolds numbers up to Re = 

10 , laminar flow was observed in the blowing section and 
was retained in the suction section for axial Reynolds 
numbers less than Re = 12000. The transition from laminar 
to fully turbulent was predicted in the entrance region of 
the suction section for axial Reynolds numbers greater than 
Re = 12000. Unlike flow in impermeable tubes, laminar 
flow was maintained for axial Reynolds numbers greater 
than 2000. For the supersonic case in the suction and 
blowing section the flow remained laminar until a shock 
wave occurred and then turbulent flow abruptly appeared, 
which showed that no transition region existed. 

Since the mathematical model is one— dimensional, 
proper expressions for the friction coefficient are necessary 
to take into account the frictional losses. Bowman (1987) 
carried out eleven numerical simulations by using the 
two-dimensional numerical model to evaluate the friction 
coefficients according to the characteristics of the vapor 
flow. The following expressions for the friction coefficients 
including the compressibility effect were correlated by using 
the friction coefficients obtained from the two-dimensional 
model. The friction coefficient for laminar flow in the 
condenser or evaporator is 

f = (1.2337 - 0.2337e _0 0363Re o)e 1 - 2M2 (14) 


The absolute value of the radial Reynolds number at the 
wall, Re 0 , is used in the evaporator and condenser. For the 

adiabatic section where the wall radial Reynolds number is 
zero, equation (14) is identical to that for the impermeable 
circular tube. For fully-developed turbulent flow in the 
condenser, the friction coefficient is 


f t = 


0.046 
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(15) 


For transition flow in the condenser entrance region, the 
friction coefficient is defined as 


f=f t + (f e -f t )e 


-0.4 12x 2 


(16) 


For the one— dimensional numerical model presented 
in this paper, the vapor flow in the evaporator and adiabatic 
sections is assumed to be laminar. In the condenser, 
laminar flow is assumed for the axial Reynolds numbers 
below Re = 12000 at the entrance of the condenser and for 
supersonic flow. When the axial Reynolds number is greater 
than Re = 12000, transition and turbulent flows are 
considered in the condenser. Also, after a shock wave 
turbulent flow is assumed in the condenser. Equation (14) 
is used to evaluate the friction coefficients in the evaporator 
and adiabatic sections and is also employed in the condenser 
when the axial Reynolds number at the entrance of the 
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condenser is less than Re = 12000 and the vapor flow is 
supersonic. When the axial Reynolds number is larger than 
Re = 12000, equation (16) is applied to the transition 
region, which is assumed to exist from the entrance of the 
condenser to about 60 percent of the condenser length based 
on the experimental data shown in Fig. 3.8 given by 
Bowman (1987). Equation (15) is used for the turbulent 
flow in the condenser. The initial location of the transition 
region, x. is equal to the location at the entrance of the 

condenser, and from Fig. 3.8 given by Bowman (1987), Xp __ 
3 y 4 and Xp _ were estimated to be 0.7 and 0.6, 
respectively. 


NUMERICAL FORMULATION 


From the many schemes (Anderson et al., 1984) 
available for the solution of the compressible flow problem, 
the Beam— Warming finite difference scheme is chosen to 
transform the governing equation (1) to the finite difference 
formulation. This scheme is a non-iterative implicit method 
and is similar to ADI for multidimensional flow problems bv 
using the factorization which retains the tridiagonal block 
matrix. 

The spatial derivatives are approximated by using 
the three-point second-order accurate central difference 
approximation for the interior points and the one-sided 
second-order accurate difference approximation for the 
boundary nodes. After the approximation operators are 
applied, the system of equations becomes 


over three grid points. The system of equations, (17), has 
the following block tri diagonal structure: 

[JKL] {tf n Dj} = {RHSj"} ( 18 ) 

where [JKL] represents the banded coefficient matrix of 
which components are 3x3 matrices for one-dimensional 

vapor flow, and {^D-} and {RHS^ n } are column vectors. 

The tridiagonal block matrix size is now (3 x I max _ 2 ) x ( 3 

x I where I is the number of nodal points. This 

max-2' max 

system of equations can be solved using the conventional 
methods for solving block tridiagonal systems of equations. 
The vector of unknowns at n + 1 time step is then 

determined by simply adding to the value of D n at n 
step. The primitive variables ( p , U, P, T) can be obtained 

from D n ^**. 

A total of 80 nodes, which is the minimum number of 
nodes to obtain accurate results, are evenly spaced in the 

axial direction and a time step of 0.1 x 10 s is used for the 
simulated heat pipe vapor flow. Since the heat flux at the 
wall in the evaporator is different from that in the 
condenser due to the different lengths for the cylindrical 
heat pipe, the coarse nodal system presented some difficulty 
to reach the steady state. For the cylindrical heat pipe, 200 
evenly-spaced nodes are used in the axial direction and a 
—3 

time step of 0.1 x 10 s is employed. 


[J W ] ^D,_i + [Kj] «“D. + (L i+1 ]« n D i+1 = [RHS]" 

' = 2 >- I max-l < 17 ) 


where 

= D n+1 - D n 


= - ttd 2 [sk ( ~ [ p l n + [ R / ) + 



where [I] is the unit matrix, [A], |P], [R] and [SJ are the 
Jacobian matrices, and e is the coefficient of the dissipative 
term to damp the oscillation. 

This time difference formula reproduces many 
different schemes with the appropriate choice of 0y and 0 

The scheme is second— order accurate in time when 0 j = 1/2 

+ ^ an( * Crst-order accurate otherwise. For Oy = 1 and 

= 1/2, the formula becomes second-order accurate in time 


RESULTS AND DISCUSSION 

1) Comparison with the Simulated Heat Pipe Vapor Flow 

A comparison of the numerical results with the 
experimental data given by Bowman (1987) is desired to 
verify the mathematical model and algorithm. However, the 
existing experimental data was obtained by simulating the 
vapor flow of a cylindrical heat pipe with a porous pipe 
which has an inside diameter of 1.65 cm and a length of 0.61 
m as shown in Fig. 1. The blowing and suction sections 


AIR IN 



| /— POROUS PIPE 

* 1 / . 

P SOURCE / 

7 — 

^ SINK 

7 

f 

V.V.V.U. JJJJJ 

.65 cm 

i 

r r r r r. a m flow^ ^ ^ 



P SOURCE 

P SINK 

— — 0.305 m — 

n | 

AIR OUT 

- — 0.305 m — - 


Fig. 1 Schematic diagram of model for air flow in the 

porous pipe. 


have equal lengths and were simulated by the injection and 
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suction of air without phase change at the interface. Also, 
uniform source and sink pressures for the blowing and 
suction regions are specified instead of the radial mass flow 
rate. Thus, to simulate the experiment equation (12) is 
used to obtain the radial mass flux with the known source 
and sink pressures. 

1.1) Transient results . The existing experimental data was 
obtained at the steady state so the present transient 
numerical results can only be compared with the transient 
numerical results for the two-dimensional model (Bowman, 
1987). For this purpose, the same geometry aaa physical 

conditions are used 6uch as the source pressure of 2.06 x 1(T* 

rt C A 

N/ar (30 psia) and the sink pressure of 1.03 x 10 N/m 
(15 psia) corresponding to case B.l. 

Initially, the velocity of the vapor is aero and the 
pressure and temperature are the same as the source 
pressure and temperature, respectively. To simulate the 
transient flow, the sink pressure is suddenly lowered to 1.03 

x 10 N/m (15 psia), while the source pressure remains the 
same as the initial pressure. This difference between the 
source and sink pressures initiates the air flow from the 
blowing section to the suction section. The pressures at the 
inlet of the blowing section, the center of the pipe, and the 
end of the suction section are plotted to compare with the 
numerical results for the two-dimensional model as shown 
in Fig. 2. 



Fig. 2 Comparison of the present numerical results 
with Bowman's 2— D model for pressure 
variations with time at three locations of the 
porous pipe: Case B.l. 


Figure 2 shows the following transient behavior of 
the air flow in the porous pipe. Since the sink pressure is 

abruptly changed from the initial pressure to 1.03 x 10 5 

N/m (15 psia) along the entire suction section, the 
pressures at the center and last nodes decrease immediately 
due to the evacuation of air. For this period, the blowing 
section pressures adjacent to the suction section start to 
decrease due to the flow of mass from the blowing section to 
the suction section, but the pressure near the beginning of 
the blowing section remains constant. Also, the mass flow 
rate from the blowing section to the suction section is not 
sufficient to influence the end of the suction section so that 
the pressure at this point decreases faster than that at the 
center of pipe. 


o 

At about 0.8 - 1.0 x 10 s, the pressure at the end 
of the suction section reaches the minimum value and then 
starts to increase while the pressure at the center of pipe 
keeps decreasing due to the frictional loss and the 
acceleration of the flow. At this time, the pressures over the 
entire blowing section become less than the initial pressure 
so that the mass flow rate is sufficient to influence the end 
of suction section. As the pressure in the blowing section 
decreases and the source pressure remains constant, the 
mass flow rate from the blowing section to the suction 
section increases. Thus, the pressure at the center node 
keeps decreasing and the pressure at the end of the suction 
section rises due to the contribution of the mass from the 
blowing section. 

—3 

At about 3.5 x 10 s, the pressures at all three 
points reach the steady state. As expected, the pressure at 
the end of the suction section does not recover completely 
due to the frictional loss at the pipe wall. Figure 2 shows 
that the present results and numerical results for the 
two-dimensional model (Bowman, 1987) are in agreement. 

1.2) Steady state results . When the present numerical 
results reach the steady state, those results are compared 
with the experimental data given by Bowman (1987). 
Numerical calculations are conducted for four different sets 
(i.e., Cases B.2, B.3, B.4, and B.5) of the source and sink 
pressures by using equations (14 — 16) for the friction 
coefficient. The cases examined are as follows: 

Case B.2: P SQ = 1.39 x 10 5 and P fik = 1.21 x 10 5 N/m 2 

Case B.3: P SQ = 2.09 x 10 5 and P $k = 1.S7 x 10 5 N/m 2 

Case B.4: P g0 = 2.68 x 10 5 and P jk = 1.69 x 10 5 N/m 2 

Case B.5: P J0 = 5.30 x 10 5 and P jk = 1.06 x 10 5 N/m 2 

Figure 3 shows the comparison of the pressure distributions 
along the axial direction. The top three lines in Fig. 3 show 



Fig. 3 Comparison of the present numerical results 
with the experimental pressure variations 
(Bowman, 1987) in the porous pipe. 


the pressure distributions for the low mass flow rates. The 
pressures in the blowing section decrease due to friction and 
the acceleration of the flow caused by mass injection, but 
the pressures in the suction section increase owing to the 
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deceleration of the flow by the extraction of mass. However, 
the pressures at the end of the suction section are less than 
those at the beginning of the blowing section because of the 
loss due to friction. In these three cases, the pressure 
distributions at the steady state correspond well to those for 
low temperature heat pipes and are in agreement with the 
experimental data. 

The fourth data profile (case B.5l shows the pressure 
variation in the axial direction for the nigh mass flow rate. 
Unlike the previous three cases, the pressure drop in the 
Wowing section is very large. The pressure ratio at the exit 
of the blowing section is about 0.4 and this ratio 
corresponds to a Mach number of M = 1. After the 
pressure decreases in the blowing section, the pressure keeps 
decreasing in the entrance region of the suction section due 
to the expansion of air even though mass removal occurs. 
Then, the pressure suddenly increases and then continues to 
increase as the flow slows down. This implies that a shock 
wave occurs at the place where the pressure changes 
abruptly. When a shock wave does not exist in the suction 
section, the pressure is supposed to decrease along the 
suction section. The one-dimensional model predicts the 
supersonic flow and shock wave in the suction section and 
the comparison of the numerical results and experimental 
data shows a good agreement except for the region 
immediately after the shock wave. 

The variations of the pressure, temperature, velocity 
and density at the steady state corresponding to cases B.4 
and B.5 are shown in Figs. 4 and 5, respectively. For the 



•s 


<• s 

“I 
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Fig. 5 Axial variations of temperature, pressure, 
velocity and density for case B.5. 


The most interesting aspect in Figs. 4 and 5 is the 
variation of temperature. The temperature at the end of the 
suction section is greater than that at the beginning of the 
blowing section. In the present model viscous dissipation is 
included. The Mach number at the entrance of suction 
section is about M=0.5 for case B.4. Since the present 
model is one— dimensional, the derivative of the axial 
velocity with respect to the radius is zero, but the friction 
effect at the interface between the wall and the vapor is 
included by using equations (14 - 16) for the friction factor. 
This effect corresponds to viscous dissipation due to the 
axial velocity derivative with respect to the radius. The 
increase in temperature is due to viscous dissipation. 

2) Comparison with Actual Vapor Flow in Cylindrical 
Heat Pipes 

The model was tested for the actual vapor flow in a 
sodium heat pipe corresponding to the experiment given by 
Ivanovskii et al. (1982) as shown in Fig. 6. The lengths of 


Q in = 560W 


Q OUT = 


Fig. 4 Axial variations of temperature, pressure, 
velocity and density for case B.4. 


low mass flow rate, the temperature and density in the 
blowing section decrease corresponding to the decrease in 
pressure and the velocity increases due to the mass injection 
as shown in Fig. 4. The Mach number, however, is less than 
M = 1 at the exit of the blowing section, so the velocity in 
the suction section decreases because of the extraction of 
mass. Also, the temperature and density increase in the 
suction section. Figure 5 shows the axial variations of the 
temperature, pressure, velocity and density for the hieh 
mass flow rate. After the sonic velocity is reached at the 
exit of the blowing section, the velocity keeps increasing in 
the suction section until a shock wave occurs. Then, the 
velocity decreases to U = 0 at the end of the pipe. As 
shown by the density profile, the air expands near the 
entrance region of the suction section and then the density 
suddenly increases after the shock wave. 



Fig. 6 Schematic diagram of model for the sodium 
vapor flow in the heat pipe. 
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the evaporator, adiabatic, and condenser sections are 0.1, 
0.05, and 0.35 m, respectively. The diameter of the vapor 
space is 0.014 m. The experimental data (Ivanovskii et aL, 
1082) represents only the saturation temperature and heat 
transfer rate (Q = 560 W) at the steady state. The initial 
conditions and actual boundary conditions applied on the 
surface of the evaporator and condenser are unknown. The 
present numerical model is transient, however, so the initial 
conditions and history of heat input and output in the 
evaporator and condenser sections are needed. 

For the numerical calculations, the uniform initial 
temperature of 810 K is used and the vapor is assumed to be 
saturated at the initial temperature. Initially, the velocity 
of the vapor is zero. The same amount of heat (Q = 560 Wj 
is uniformly applied on the surface of the evaporator and the 
convective boundary condition is used on the surface of the 
condenser. The reference temperature of 300 K for the 
convective boundary condition is employed and the heat 
transfer coefficient is determined by iteration. At first, an 
arbitrary initial heat transfer coefficient is guessed. When 
the numerical results reach the steady state, the saturation 
temperature at the end cap of the evaporator is compared 
with the experimental data at the same location. This 
procedure is repeated until the same temperature is 
obtained at the steady state. For this test, the heat transfer 
2 

coefficient of 69.1 W/m — K is used. 

Figure 7 shows the axial variation of the saturation 
and vapor temperatures, pressure, velocity and density 
obtained from the present numerical model and the 
experimentally measured saturation temperature 
distribution (Ivanovskii et al., 1982). The pressure, 



Fig. 7 Axial variations of temperature, pressure, 
velocity and density of the sodium heat pipe 
at steady state. 


temperature and density in the evaporator decrease and the 
velocity increases due to the injection of mass and the effect 
of friction. In the adiabatic section, the vapor temperature 
increases because of viscous dissipation and the pressure 
decreases owing to friction at the interface. The density 
also decreases while the velocity continues to increase. The 
vapor temperature in the condenser continues to increase 
due to viscous dissipation so that the vapor temperature at 
the end cap of the evaporator is less than that at the end 
cap of the condenser. The pressure recovery in the 
condenser is almost negligible. This may result from 
dominant friction effect at the interface of the condenser 
compared to the effect of mass extraction in this long 


condenser. In the adiabatic section a difference between the 
calculated and measured saturation temperatures is 
observed, but the trend of the saturation temperature 
variation is the same. The trend of the vapor temperature 
in the condenser is quite different from that of the 
saturation temperature in the same region, so the saturation 
temperature may not be assumed to be the vapor 
temperature in the condenser for the one-dimensional 
model. 

The effect of viscous dissipation in the vapor flow for 
the high temperature heat pipe is investigated as shown in 
Fig. 8. The vapor temperature at the end cap of the 
evaporator with viscous dissipation is less than that at the 
end cap of the condenser. However, the vapor temperature 
at the end cap of the evaporator without viscous dissipation 
is almost the same as that at the end cap of the condenser. 
In the adiabatic section, the vapor temperature with viscous 
dissipation increases while that without viscous dissipation 
does not change. Therefore, the viscous dissipation should 
be taken into account in the energy equation for the high 
temperature heat pipe. 


d 



Heat Pipe Axial Distance, m 


Fig. 8 Comparison of the temperature variations 

with and without viscous dissipation for 
actual vapor flow. 


CONCLUSIONS 

A model for the transient one-diraensional 
compressible vapor flow in the cylindrical heat pipe is 
developed. This model predicts the vapor flow in cylindrical 
heat pipes as well as simulated heat pipes for the subsonic, 
sonic, and supersonic flows under transient and steady state 
conditions. The vapor flow quickly reaches the steady state 
condition. The distributions of the temperature and 
pressure during the transient state are quite different from 
those for the steady state. The viscous dissipation terms 
play an important role in the energy equation and have to 
be taken into account. The one-dimensional compressible 
model predicts the experimental data well for the cylindrical 
heat pipe and the Simulated heat pipe at the steady state. 
The experimental data for the transient state arc needed to 
understand clearly the transient behavior of the vapor flow 
in the heat pipe both at low and high temperatures. The 
one-dimensional model can reduce the computational effort 
needed to solve the vapor flow problem. 
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